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Abstract  -  Electroencephalogram  (EEG)  analysis  remains 
problematic  due  to  both  lack  of  understanding  of  the 
origins  of  the  signal  and  inadequate  evaluation  methods. 
In  spite  of  these  shortcomings,  the  EEG  is  a  valuable  tool 
in  the  evaluation  of  some  neurological  disorders  as  well  as 
in  the  evaluation  of  overall  cerebral  activity.  It  becomes 
more  useful  when  combined  with  other  clinical 
parameters.  The  focus  of  the  work  described  here  is  two¬ 
fold.  New  chaotic  methods  are  introduced  for  EEG 
evaluation  coupled  with  a  hybrid  system  approach  that 
permits  the  combination  of  the  EEG  results  with  clinical 
parameters  to  form  a  comprehensive  decision  model.  The 
system  is  illustrated  in  an  application  for  diagnosis  of 
dementia.  Extensions  can  easily  be  made  to  applications 
such  as  evaluation  of  brain  activity  during  surgery. 
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I.  Introduction 

The  exact  mechanism  of  the  generation  of 
electroencephalogram  (EEG)  signals  is  not  understood,  due  in 
part  to  the  lack  of  appropriate  theoretical  models  and 
appropriate  measurements  to  adequately  describe  and  dissect 
the  EEG  signals.  Basic  approaches  to  signal  analysis  have 
relied  on  Fourier  analysis,  cross-correlation,  auto  correlation, 
and  other  techniques  to  determine  if  the  signal  is  stationary 
[1],  While  these  approaches  have  proved  useful  in  many 
areas,  analysis  of  many  medical  time  series  such  as  ECGs  and 
EEGs  are  still  problematic.  Conventional  EEG  evaluation 
methodologies  are  useful  but  limited  and  potentially 
problematic  form  both  theoretical  and  practical  standpoints. 
EEG  signals  are  considered  as  the  results  of  the  combined 
dynamic  activity  of  neuronal  populations.  Models  including 
excitatory  and  inhibitory  circuits  with  feedback  loops  have 
been  adopted  to  explain  the  oscillation  property  of  EEG 
activity  [2],  Clinical  correlations  of  the  dominant  signal 
frequencies  and  visual  detection  of  paroxysmal  events  such  as 
spikes  or  sharp  waves  have  been  the  mainstay  of  clinical 
neurological  interpretation  of  EEG  recording.  The  traditional 
approach  to  EEG  analysis,  Fourier  analysis  provides  a 
quantitative  tool  to  examine  signal  frequencies  and  relative 
loads.  It  is  almost  certain  that  conventional  Fourier  analysis 
cannot  represent  the  entire  spectrum  of  biological  activities. 
In  addition,  some  of  the  assumptions  such  as  the  stationarity 
of  the  signal  are  not  valid.  Signal  averaging  and  analysis 
based  on  short  intervals  ranging  from  one  to  four  sections  are 
inadequate.  These  problems  may  bias  the  analysis  [3], 


Clinical  utility  of  the  EEG  is  also  limited  by  the  frequent 
lack  of  specificity  of  the  EEG  abnormality.  Generalized 
slowing  during  an  EEG  tracing  unrelated  to  drowsiness  can 
be  an  indication  of  generalized  cerebral  dysfunction  due  to 
metabolic  derangement,  neurodegenerative  disorders,  or 
infectious  or  inflammatory  diseases.  Conventional  EEGs 
include  18  channels  with  only  limited  resolution  for 
localization,  imposing  yet  another  limitation.  More 
comprehensive  linear  and  nonlinear  analyses  of  the  EEG 
signals  described  here  not  only  have  practical  utility  [4]  but 
can  also  open  new  windows  for  studying  the  significance  of 
the  EEG  signal  in  the  understanding  of  the  basic 
neurophysiological  functioning  of  the  human  cerebral  cortex. 
A  nonlinear  approach  using  continuous  chaotic  modeling  that 
provides  measurements  of  the  level  of  variability  of  the  EEG 
is  described  below.  The  method  is  illustrated  in  an 
application  for  diagnosis  of  dementia  and  can  be  extended  to 
analysis  of  brain  activity  during  surgery. 

II.  Methodology 

A.  Theoretical  Basis  for  Chaotic  Analysis  of  Time  Series 

The  basic  common  thread  in  chaos  theory  is  the  recursive 
evaluation  of  seemingly  simple  functions  that  produce 
unexpectedly  complex  results.  An  iterative  function  does  not 
suddenly  become  chaotic,  but  rather  goes  from  the  stage  of 
convergence  to  a  single  value  to  a  bifurcation,  or  convergence 
to  two  values.  Additional  bifurcations  occur,  and  finally 
chaos  results.  As  an  example,  consider  logistic  equation 

an=  A  an_j(  1  -  an_i)  2<A<4  (1) 

where  A  is  a  constant  whose  value  changes  the  behavior  of 
the  function.  The  recursion  is  dependent  on  the  selection  of 
a0,  which  must  be  chosen  between  0  and  1.  For  increasing 
values  of  A,  the  equation  progresses  from  single  value 
convergence  to  chaos.  Within  the  chaotic  area,  regions  of 
stability  unexpectedly  appear.  For  integer  values  of  n,  this 
function  exhibits  chaotic  properties  for  A  >  3.57.  These 
properties  include  apparent  lack  of  periodicity  and  sensitivity 
to  initial  conditions.  The  picture  changes,  however,  if 
continuous  values  instead  of  integers  are  considers.  The 
exact  solution  of  (1)  at  A  =  4  is: 

an  =  i/2[l-T2n(l-2a0)]  (2) 

where  Tn(x)  is  the  Chebyshev  function  [5],  The  authors  have 
derived  a  soft  solution  for  (1)  for  all  values  of  A  [6]: 


Report  Documentation  Page 


Report  Date  Report  Type  Dates  Covered  (from...  to) 

25  Oct  2001  N/A 

Title  and  Subtitle  Contract  Number 

EEG  Analysis  Based  on  Chaotic  Evaluation  of  Variability 

Grant  Number 
Program  Element  Number 

Author(s)  Project  Number 

Task  Number 
Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es)  Performing  Organization  Report  Number 

California  State  University  Fresno,  CA  93740 

Sponsoring/Monitoring  Agency  Name(s)  and  Address(es)  Sponsor/Monitor’s  Acronym(s) 

US  Army  Research,  Development  &  Standardization  Group  - 

(UK)  PSC  802  Box  15  FPO  AE  09499-1500  Sponsor/Monitor’s  Report  Number(s) 

Distribution/ Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

Papers  from  23rd  Annual  International  Conference  of  the  IEEE  Engineering  in  Medicine  and  Biology  Society,  October 
25-28,  2001,  held  in  Istanbul,  Turkey.  See  also  ADM001351  for  entire  conference  on  cd-rom. 

Abstract 


Subject  Terms 


Number  of  Pages 

4 


Assume  a  solution  of  the  type 


/ 

an=  £  akTk(2nx) 
k=0 


(3) 


where  Tk(x)  is  the  Chebyshev  function  of  the  first  kind  and  n 
is  a  real  number.  We  assume  /  to  be  the  number  of  points  in 
the  interval  0  <  n  <  1 .  Thus 


/ 

an2  =  £  ak2  Tk2  (2“x)  +  2  £  a;  aj  T^"  x)Tj  (2"  x)  (4) 

k=0  j  >  i 

i=0. 1 _ _  /-I 

1=1,2,...,  / 


The  conjecture  adopted  is  that  going  from  one  point  to 
another  implies  adding  a  Chebyshev  polynomial.  Hence 

21 

an+i  =  £  bkTk(2n  x)  (5) 

k=0 


where  n  is  assumed  to  be  a  real  number.  By  imposing 
appropriate  boundary  conditions  one  obtains  a  unique 
solution  to  these  nonlinear  equations  involving  300  variables. 
Values  for  n  >  1  are  obtained  by  applying  the  logistic 
equation  to  the  points  obtained  for  0  <  n  <  1 . 

Plots  of  the  continuous  solution  show  no  dramatic  change 
in  behavior  at  A  =  3.57,  but  rather  a  well-defined  increase  in 
variability,  as  illustrated  by  examining  the  second-order 
difference  plots  generated  by  the  soft  solution.  Second  order 
difference  plots  are  generated  by  plotting  an+2  -  an+i  versus 
an+i  -  an,  where  an  is  the  value  of  the  time  series  at  time  n. 
Theoretical  plots  at  A  =  3.57  and  A  =  4.0  are  shown  in  Fig.  1. 
Equivalent  plots  can  be  used  for  the  evaluation  of  time  series, 
as  shown  in  Fig.  2  for  electrocardiogram  (ECG)  evaluation  of 
a  normal  patient  and  a  patient  with  congestive  heart  failure. 
In  this  approach,  rather  than  defining  a  time  series  as  chaotic 
or  not  chaotic,  it  is  evaluated  in  terms  of  degree  of  variability 
or  chaos.  To  quantify  the  level  of  variability,  the  central 
tendency  measure  (CTM)  is  used,  which  is  computed  by 
selecting  a  circular  region  around  the  origin  of  radius  r, 
counting  the  number  of  points  within  the  radius,  and  dividing 
by  the  total  number  of  points  t.  Then 

t-2 

n=  [£  8  (dj)]/(t-2)  (6) 

i=l 

where  8(d;)  =  1  if  [(ai+2-ai+1)2+(ai+i-ai) 2]'5  <  r  and  0  otherwise. 

The  CTM  measure  has  been  shown  to  be  effective  in 
analysis  of  ECG  data  in  several  ways:  as  an  independent 
measure,  combined  with  other  ECG  measures  in  a  neural 
network  model,  and  combined  with  clinical  parameters  in  a 
neural  network  model  [7],  Preliminary  studies  have  shown 
its  feasibility  for  use  in  EEG  analysis  [8], 
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Figure  1 :  Theoretical  Second-Order  Difference  Plots 
(a)  A  =  3.57  and  (b)  A  =  4.0 
x  —  an+i  -  an,  y  =  an+2  -  &n+i 


(a)  (b) 

Figure  2:  Second-Order  Difference  Plots  for  a  Normal  Subject  (a)  and  a 
Congestive  Heart  Failure  Subject  (b) 


B.  Implementation  for  Rapid  Evaluation 

Another  problem  encountered  with  EEG  analysis  is  the 
large  number  of  data  points.  For  a  10-minute  run  for  each 
channel,  approximately  75,000  points  are  recorded,  with  a 
minimum  of  eighteen  channels.  In  an  application  such  as 
surgery,  any  useful  evaluation  must  be  capable  of  running  in 
real  time.  While  plotting  a  second-order  difference  plot  is 
time-consuming,  the  CTM  measure  can  be  done  directly  and 
rapidly  enough  to  produce  real-time  results 

C.  Preprocessing  for  EEG  Signal  Analysis 

EEG  signals  require  preprocessing  for  removal  of  noise 
and  identification  of  peaks.  Two  methods  that  are  currently 
being  tested  are  the  use  of  a  peak  identification  algorithm 
developed  by  the  authors  [9]  and  wavelet  processing  that 
allows  the  identification  of  peaks  of  varying  amplitudes  [10]. 

D.  Hybrid  System  for  Data  Analysis 

The  hybrid  system  Hypermerge  is  used  to  combine 
multiple  EEG  results  with  clinical  parameters  [11], 
Hypermerge  has  three  components: 

Rule -based  component  (EMERGE) 

Data-based  component  (Hypernet) 

Chaotic  Analysis  of  Time  Series  (CATS). 


In  this  application,  the  rule-based  component  is  used  to 
include  expert  opinion,  the  neural  network  model  is  used  to 
combine  EEG  summary  results  with  clinical  parameters  and 
neuropsychological  testing  results,  and  the  chaotic  analysis  is 
used  for  evaluation  of  the  EEG.  The  knowledge-based 
component  can  include  a  wide  range  of  information,  ranging 
from  impressions  of  mental  status  to  human  interpretation  of 
imaging  and  EEG  results. 

III.  Results 

A.  Collection  ofCTM  Variability  Data 

Preliminary  EEG  data  has  been  collected  to  evaluate  the 
feasibility  of  this  approach.  Data  were  collected  at  a  rate  of 
250  samples/second  with  a  periodic  2-second  delay  for 
storage  requirements.  Digital  EEG  runs  lasted  approximately 
10  minutes  and  consisted  of  approximately  75,000  points. 
Each  data  point  consists  of  a  consecutive  number  and  two 
channels  of  output.  The  output  value  for  each  channel  is  a 
positive  or  negative  integer  indicating  the  current  amplitude. 
Two  channels  are  selected  from  the  21  available  for  this 
preliminary  analysis.  The  channels  selected  for  recording 
were  T3-T5  and  T4-T6.  T3  to  T6  are  based  on  standard  EEG 
electrode  placement.  T3-T5  locates  over  the  left  temporal 
area  with  T4-T6  over  the  right. 

B.  Comparison  of  Lead  Activity 

Symmetry  is  always  an  important  medical  indicator  of 
abnormal  states.  This  is  especially  true  in  symmetric  organs 
such  as  the  brain  and  is  often  used  in  the  evaluation  of 
medical  images.  It  can  also  be  used  in  the  evaluation  of 
EEGs.  Electrodes  are  placed  to  cover  each  of  the  different 
lobes  of  the  brain,  frontal,  parietal,  occipital,  and  temporal, 
and  are  symmetrically  placed  on  the  right  and  left  sides,  as 
described  in  the  data  collection  above.  Thus  the  summary 
measures  for  corresponding  areas  can  be  compared  to  each 
other  to  see  if  a  similar  level  of  activity  is  occurring. 

Table  I  shows  CTM  measures  based  on  two  methods  of 
analysis  for  leads  placed  on  the  left  and  right  temporal  lobes 
for  6  Alzheimer’s  patients  and  2  normal  controls.  Two 
methods  of  analysis  were  used: 

The  second-order  difference  plot  generated  based  on  each 
point  in  the  time  series:  This  analysis  is  based  on 
amplitude  values  indicating  the  level  of  electrical  activity. 
Time  between  peak  occurrence:  Time  between  peaks  was 
used  as  an,  the  nth  point  in  the  series,  to  generate  the 
second-order  difference  plot.  This  analysis  is  based  on 
frequency  values  of  the  occurrence  of  the  peaks  as 
determined  by  the  peak  identification  algorithm. 

C.  Hybrid  Evaluation 

Hypermerge,  the  hybrid  system  developed  by  the  authors, 
is  used  to  implement  the  decision  strategy.  The  knowledge- 
based  component  and  neural  network  are  described  briefly. 


TABLE  I 

COMPARISON  OF  SYMMETRIC  LEADS 


ID# 

Left  Temporal 

Am.  (r=0.1)  Fr.  (r=05) 

Right  Temporal 

Am.  (r=0.1)  Fr.  (r=.05) 

N, 

0.54 

0.29 

0.50 

0.31 

n2 

0.59 

0.57 

0.59 

0.55 

Ai 

0.81 

0.52 

0.80 

0.53 

A2 

0.40 

0.62 

0.41 

0.62 

a3 

0.67 

0.28 

0.66 

0.26 

A4 

0.60 

0.44 

0.60 

0.43 

A5 

0.68 

0.40 

0.68 

0.40 

A6 

0.58 

0.18 

0.58 

0.22 

Nj:  Normal  controls;  A 

i!  Alzheimer’ s  Patients 

Am.: 

Amplitude;  Fr.:  Frequency 

Knowledge-Based  Component  ( Emerge ) 

The  knowledge  base  is  in  the  form  of  rules  that  are  derived 
from  expert  input.  An  example  of  a  rule  used  to  supplement 
EEG  analysis  for  the  diagnosis  of  dementia  is  given  below: 


Rule  # 

Premise 

Weight 

IF 

EEG  shows  low  variability  in  temporal  lobe 

0.4 

Score  <  1 1  on  MMSE 

0.4 

Family  hx  of  Alzheimer’s 

0.2 

THEN 

Prescribe  aricept 

Threshold  0.6 

The  rule  is  evaluated  using  approximate  reasoning 
techniques  to  aggregate  evidence  [12],  In  addition  to  the 
weighting  factors,  the  user  determines  the  level  of  presence  of 
each  condition  by  entering  a  value  between  zero  and  ten 
indicating  the  severity.  The  value  is  then  divided  by  ten, 
yielding  values  between  0  and  1,  inclusive.  The  weighting 
factors  on  the  premises  are  normalized  to  sum  to  1.  If  the 
objective  is  to  evaluate  brain  function  during  surgery,  a  new 
rule  base  is  substituted.  An  important  feature  of  this  design  is 
that  no  algorithmic  changes  need  to  occur  when  the  new  rule 
base  is  inserted.  Evaluation  of  the  rules  for  the  specific 
patient  case  is  very  rapid  and  can  easily  be  done  within  the 
real  time  requirement  for  surgery. 

Neural  Network  Model  (Hypernet) 

The  neural  network  model  is  used  in  two  ways:  to 
combine  different  measures  of  variability  obtained  from  the 
chaotic  analysis  and  to  combine  the  chaotic  measures  with 
other  clinical  and  neuropsychological  parameters.  Variables 
for  input  nodes  along  with  their  sources  for  dementia 
evaluation  are  given  in  Table  II.  If  the  application  is  adjusted 
for  surgery,  the  first  12  parameters  are  still  relevant.  Nodes 
13-20  may  also  be  relevant.  The  clinical  inputs  can  be 
replaced  with  any  pertinent  clinical  components.  A  new 
decision  model  is  then  easily  obtained  using  the  Hypernet 
learning  algorithm.  Once  the  decision  model  is  established, 
computation  of  patient-specific  input  values  can  easily  be 
accomplished  in  real  time. 


Table  ii 

INPUT  TO  NEURAL  NETWORK  MODEL 


Acknowledgment 


Node  #s 

Contents 

Origin 

NrN8 

CTM 

Chaotic  Analysis 

N9-N12 

CTM  Difference 

Chaotic  Analysis 

N13-N20 

Activity  levels 

Functional  Imaging 

N2i 

MMSE 

Neuropsychological  Test 

N22 

Genetic  Factor  (y/n) 

Genetic  Testing 

N23 

Family  Hx  (y/n) 

Interview 

N24 

Visible  impairment  (y/n) 

Exam 

Location  Codes 

T  (temporal),  P  (parietal).  O  (occipital).  F  (frontal) 

L  (left),  R  (right) 

Node  Definitions 

NrN8 

TL,  TR,  PL,  PR,  OL,  OR,  FL,  FR 

N9-Ni2 

T,  P, O,  F 

N13-N20 

TL,  TR,  PL,  PR,  OL,  OR,  FL,  FR 

All  variables  are  continuous  unless  otherwise  indicated. 


IV.  Discussion 

The  approach  to  continuous  chaotic  modeling  described 
here  has  previously  been  shown  to  be  useful  in  differentiation 
of  categories  of  cardiac  disorders  using  ECG  analysis.  In 
these  applications,  the  combination  of  the  chaotic  analysis 
with  clinical  parameters  through  the  use  of  a  neural  network 
decision  model  increased  sensitivity,  specificity,  and 
accuracy  of  results.  In  preliminary  results  in  diagnosis  of 
dementia,  the  same  approaches  look  promising  for  the 
analysis  of  EEGs.  The  same  paradigm  is  followed  in  which 
the  chaotic  analysis  is  used  as  part  of  a  more  general  neural 
network  model.  As  a  third  step,  these  two  modalities  are 
included  in  a  comprehensive  hybrid  system  that  also  permits 
the  inclusion  of  expert-derived  knowledge  through  the  use  of 
rules.  The  intent  of  the  initial  implementation  that  includes 
EEG  analysis  is  to  diagnose  dementia  at  an  early  stage  at 
which  the  potential  for  successful  treatment  is  greater. 
However,  the  model  is  easily  adaptable  to  other  applications 
of  EEG  analysis,  including  monitoring  during  surgery.  The 
patient-analytic  parts  of  the  system  have  been  implemented  to 
run  in  real  time  so  that  use  during  surgery  is  feasible.  The 
EEG  analysis  for  a  surgery  application  remains  the  same. 
The  neural  network  model  only  requires  retraining  on  a  new 
data  set  after  appropriate  inputs  have  been  determined.  The 
knowledge  base  component  requires  simple  replacement  of 
the  rule  base.  None  of  the  three  components  of  the  hybrid 
system  require  any  alteration  for  the  new  application. 

V.  Conclusion 


This  work  was  supported  by  NIH  Grant  MBRS -SCORE 
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Analysis  using  chaotic  parameters  presents  a  novel 
approach  for  the  extraction  of  information  from  the  complex 
mix  of  signals  that  make  up  the  EEG.  Early  work  using  both 
frequency  and  amplitude  analysis  looks  promising.  Use  of 
the  chaotic  parameters  in  a  comprehensive  decision  model 
expands  the  potential  for  using  the  EEG  as  an  important 
clinical  parameter  in  the  diagnosis  of  disease  and  in  the 
monitoring  of  patients. 


